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ABSTRACT 

We combine new sensitive, wide-field CO data from the HERACLES survey with ultraviolet and 
infrared data from GALEX and Spitzer to compare the surface densities of H2 , Sh2 , and the recent 
star formation rate, Ssfr, over many thousands of positions in 30 nearby disk galaxies. We more 
than quadruple the size of the galaxy sample compared to previous work and include targets with 
a wide range of galaxy properties. Even though the disk galaxies in this study span a wide range 
of properties, we find a strong, and approximately linear correlation between Ssfr and I]h2 at our 
common resolution of 1 kpc. This implies a roughly constant median H2 consumption time, r^^p = 
Eh2/5]sfr, of ~ 2.35 Gyr (including heavy elements) across our sample. At 1 kpc resolution, there is 
only a weak correlation between Sh2 and Tp^p over the range Eh2 ~ 5-100 Mq pc~^, which is probed 
by our data. We compile a broad set of literature measurements that have been obtained using a 
variety of star formation tracers, sampling schemes and physical scales and show that overall, these 
data yield almost exactly the same results, although with more scatter. We interpret these results 
as strong, albeit indirect evidence that star formation proceeds in a uniform way in giant molecular 
clouds in the disks of spiral galaxies. 

Subject headings: galaxies: evolution — galaxies: ISM — radio lines: galaxies — stars: formation 



1. INTRODUCTION 

Giant molecular clouds (GMC s) are the sites of star 
formation in the Milky Way fe.g.. lBlitd[T993l ). Therefore, 
it is not surprising that a strong correlation is observed 
between tra cers of molecular gas and recent star for- 
; mation (e.g ,jRownd & Young"1999'; ' Wong fc Blitg[2002l: 
. ILerov et all [2008; Bigicl ct al. 2008), while the corre- 
lation between atomic gas and recent star formation 
is found to be weak or absent within galaxies (e.g., 
. iKennicutt eFall [20071: iBigiel et al.ll2Q08l ). The details of 
this correlation have important implications. Its evo- 
lution over cosm ic time informs our understanding o f 
I galaxy assembly (jDaddi et al.ll2010l : iGenzel et aLll2010f ). 

The finding of short molecular gas consumption times 
I compared to galaxy lifetimes highlights the importance of 
fueling the inner disks of galaxies. The relatively low effi- 
ciency of star formation per dynamical time requires that 
the star formation process be more complex than simple 
gravitational collapse (e.g., McKe_e_& Ostrikcr 2007). Fi- 
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nally, the relationship between star formation and molec- 
ular gas is an important input and benchmark for models 
attempting to reproduce today's galaxies or galaxy pop- 
ulations. 

The importance of this topic has led to several stud- 
ies of the relationship between surface densities of H2 
and the star format ion rate. Many of these focus on 
single galaxies (e.g., 'Hever et al.l 1200 4": 'Kennicutt et al. 
2007; Schuster ct al. 2007; Bl anc et al . 2009; Verlcv et al. 
2010; Rahman ct al. 2010) or a small sample (e.g. , 
Wong fc Blitz 2002: Lerov ct al. 2008: Bigiel et aT1[2008l : 
Wilson et al.l 120091: IWarren et al.l l2010r ). Restricted 
by the availability of sensitive and wide-field molecu- 
lar gas maps, stu d ies of larg e sets of galaxies (e.g., 
lYoung et ajl Il996t 'Kennicutt '1998; 'Rown d fc Youngl 
1999; MurgiaeLai]l2002; Lcroy ct al. 2005) mostly used 
integrated measurements or a few pointings per galaxy. 
To date, no homogeneous analysis of the correlation be- 
tween the star formation rate and H2 surface densities in 
a large set of nearby galaxies at good spatial resolution 
exists. 

In this letter we take this next logical step, compar- 
ing molecular gas — traced by CO emission — to re- 
cent star formation — traced by ultraviolet and infrared 
emission — at 1 kpc resolution across a large sample of 
30 nearby galaxies. This samp le is signifi c antly larger 
and more diverse than that of Bigie l et al.l ([20081 here- 
after B08). From 2007-2010, the H:ERA CO-Line Extra- 
ga lactic Survey (HERACLES, first maps are presented 
in ILerov et"ani2009) collaboration used the IRAM 30-m 
telescop^m to construct maps of CO J = 2 ^ 1 emis- 
sion from 48 nearby galaxies. Because the target s over- 
lap surveys by Spitzer (mostly SINGS, Kennicut t et al.l 

" IRAM is supported by CNRS/INSU (France), the MPG (Ger- 
many) and the IGN (Spain). 
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TABLE 1 



Galaxy Sample 


Galaxy 


D 
[Mpc] 


Galaxy 


D 
[Mpc] 


NGC OSST^ 


24.7 


NGC 4254^^ 


20.0 


NGC 0628^°* 


7.3 


NGC 4321 


14.3 


NGC 0925 


9.2 


NGC 4536 


14.5 


NGC 2403 


3.2 


NGC 4559 


7.0 


NGC 2S41 


14.1 


NGC 4569"^ 


20.0 


NGC 2903 


8.9 


NGC 4579^^ 


20.6 


NGC 2976 


3.6 


NGC 4625 


9.5 


NGC 3049 


8.9 


NGC 4725 


9.3 


NGC 3184^08 


11.1 


NGC 47362"** 


4.7 


NGC 3198 


13.8 


NGC 5055'^''** 


10.1 


NGC 3351 


10.1 


NGC 5194^08 


8.0 


NGC3521B08 


10.7 


NGC 5457 


7.4 


NGC 3627 


9.3 


NGC 5713^1 


26.5 


NGC 3938 


12.2 


NGC 6946'3"8 


5.9 


NGC 4214 


2.9 


NGC 7331 


14.7 



^'^^ Target from B08. 

Too distant to reach 1 kpc resolution, included 
in the 1 kpc plots at their respective native reso- 
lution. 



20M) and GALEX (mostly the NGS, iGil de Paz eTall 
2007), excellent multiwavelength data are available for 
most targets. 

2. METHOD 

We study all galaxies meeting the following criteria: 1) 
a HERACLES map containing a robust CO J — 2 1 
detection, 2) GALEX far UV (FUV) and Spitzer in- 
frared data at 24/im (IR), and 3) an inclination < 75°. 
The first condition excludes low mass galaxies without 
CO detections. The second removes a few targets with 
poor Spitzer 24^m data. The third disqualifies a hand- 
ful of edge-on galaxies. We are left with 30 disk galax- 
ies, listed in Table P along with distances adopted from 
iWalter et all (120081 ^ LEDA, and NED. This sample is 
more than four times larger than that of BOS and spans 
a substantial range in metallicities (8.36 < z < 8.93^3 
and mass (8.9 < \og{A'Q < 11.0^3- 

We follow the approach of B08 with only a few modifi- 
cations. B08 compared the first seven HERACLES maps 
to FUV, IR, and Ha emission to infer the relationship be- 
tween the surface density of H2 , Eh2 , and the star forma- 
tion rate surface density, Esfr- As in B08, we estimate 
SH2 from HERACLES CO J = 2 ^ 1 emission. We 
assume a Galactic Xco — 2 x 10^° cm^^ (K km s^^)^^, 
correct for inclination, include helium in our quoted sur- 
face densities (a factor of 1.36, a difference from B08), 
and adopt a CO line ratio I {2 - 1)//(1 - 0) = 0.7 (note 
that B08 used a ratio of 0.8). 

We estimate Ssfr (inclination corrected) using a com- 
bination of FUV emission and 24/zm emission. FUV 
emission traces mainly photospheric emission from O 
and B stars, with a typical age of ^ 20-30 Myr 
(jLeitherer et al.l 119991 : iSahm et alJ l2007| ) but sensitive 
to populations up to 100 Myr of age. Infrared emis- 

Metallicities are adopted from Moust akas et all l|201(3) where 
available and supplemente d by data fr om the compilations in 
ICalzetti et all 1 1201^1) and lMarble et al.1 (I20TO|'I . 

Ste llar masses are estimat ed using the near IR luminosi- 
ties from Dale et al.l 1)20071 . [20091 1 and the mass-to-light ratio from 
ILerov et al. (,2003). 



sion at 24/zm comes from dust mainly heated by young 
stars. This emission correlates closely with other sig- 
natures of recent star formation, especially Ha emis- 
sion, and so has been used to cor rect optical and UV 
tracers for the e f fects of extinction (IC alzetti et al."2007|: 
iKennicutt et al.l l2007| ). ILerov et aLI (,200&) motivated 
this FUV-IR combination, showing that it reproduces 
other estimates of Ssfr with ^ 50% accuracy down to 
SsFR « 10-3 Mq yr-i kpc-2. 

For 24 galaxies, we use FUV maps fr om the Nearby 
Galaxy Survey (NGS, IGil de Paz eTaH Eoof). for five 
targets from the All-sky Imaging Survey (AIS) and for 
one galaxy we use a map from the Medium Imaging Sur- 
vey (MIS). We use maps of IR emission at 24/im from 
the Spitzer Infrared Nearby Galaxies Survey (SINGS, 
IKennicutt et al.ll2003l) and t he Local Volume Legacy Sur- 
vev (LVL. iDale et al.l [20091 ). Handling of the maps fol- 
lows B08. 

We convolve the IR and FUV maps to the 13" 
(FWHM) resolution of the HERACLES data. Given the 
wide distance range of our sample, 13" resolution corre- 
sponds to physical scales from 180 pc to 1.7 kpc. To avoid 
being influenced by physical resolution, we create a sec- 
ond set of maps at a common physical resolution of 1 kpc 
(FWHM), appropriate to carry out a uniform analysis. 
Five galaxies are too distant to reach 1 kpc resolution. 
We include them in our "kpc" analysis at their native 
resolution, 1.4 kpc on average (excluding them does not 
change our conclusions). 

The HERACLE S maps are rnasked to include only sig- 
nificant emission ()Lerov et al.ll2009D . The exact com- 
pleteness of each map in mass surface density depends 
on the inclination and, for fixed spatial resolution, the 
distance of the target. A typical noise level is 25 mK per 
5.2 km s^^ channel at 13" resolution. For the most dis- 
tant, face-on systems this limitis Ico > 0.8 K km or 
5 Mq pc"^ for our adopted Xco and line ratio. Closer 
or more inclined systems will be complete to lower SH2- 

We sample both sets of maps, one at 13" and one at 
1 kpc resolution, using a hexagonal grid spaced by one 
half- resolution element. We keep only sampling points 
inside the B-band 25"^ magnitude isophotal radius, r25, 
and where the HERACLES mask includes emission. At 
13" resolution, this yields Ssfr and Sh2 estimates for 
a total of ~ 27, 000 points 5, 000 independent mea- 
surements) in 30 nearby star-forming galaxies. At 1 kpc 
resolution, this number drops to ~ 12,000 (~ 2000 inde- 
pendent) measurements. 

3. RESULTS 

Figuredlshows our data in SgFR-SH2 space. The upper 
panels present measurements at a common angular reso- 
lution of 13", the lower panels show results for a common 
physical scale of 1 kpc. The left panels show contours in- 
dicating the density of data with each galaxy weighted 
equally. The right panels directly show each data point. 
Dotted lines in each plot indicate constant molecular gas 
depletion times, t^^^ = I]h2/SsfRj i-G-, fixed ratios of 
H2-to-SFR. 

To make the contour plots, we divide the I]sfr-Sh2 
space into 0.1 dex- wide cells to grid the data. During 
gridding, we assign each data point a weight inversely 
proportional to the number of data points for the galaxy 
that it was drawn from. This assigns the same total 
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Fig. 1. — Star formation rate surface density, SgpR, estimated from FUV+24/im emission as a function of molecular gas surface density, 
Sh2> estimated from CO J = 2 — > 1 emission for 30 nearby disk galaxies. The left panels show data density with equal weight given to 
each galaxy. Purple, red, orange, and green contours encompass the densest 25, 50, 75, and 90% of the data. The right panels show each 
measurement individually as a black dot. The red points indicate running medians in SsFR ^ function of Sh2 and the error bars show 
the la log-scatter in each Eh2 bin. In both panels, dotted lines indicate fixed H2 depletion times in yr. Measurements in the top panels 
are on a common angular scale of 13", those in the bottom panels are on a common physical scale of 1 kpc. All panels show a strong 
correlation between Ssfr and Sh2 with the majority of data having T^^p ~ 2.3 Gyr. 



weight to each galaxy, ensuring that a few large galaxies 
do not drive the overall distribution. Contours indicate 
the density of sampling points in each cell. 

The scatter plots on the right treat all measurements 
equally, which leads large galaxies to dominate the dis- 
tribution. While the contour plots treat a galaxy as the 
fundamental unit, the scatter plots treat each region as 
a key independent measurement. The red points show a 
running median in Ssfr as a function of Sh2- Though 
treating Eh2 as an independent variable is not rigorous, 
this binning is a useful way to guide the eye. We only bin 
where Sh2 > 5 Mq pc~^ and we are confident of being 
complete. 

All four plots reveal a strong correlation between Ssfr 
and • In this letter we focus our quantitative analysis 



on the right hand plots, which weight every measurement 
equally. The Spearman rank correlation coefficient across 
all data is r = 0.8 at 1 kpc resolution, indicating a strong 
correlation between Ssfr and Sh2- We find a median H2 
depiction time t^^^ ~ 2.35 Gyr with la scatter 0.24 dex 
(« 75%). The results at fixed 13" resolution are similar, 
median is ~ 2.37 Gyr and r = 0.7. 

It is common to parameterize relationships between 
gas and star formation using power law fits. This can be 
problematic physically, because data from widely vary- 
ing environments are often not well- described by a single 
power law (BOS. lBigiel et al.ll2010B . It is also challeng- 
ing practically, becau se of, e.g., issues of completeness 
and upper limits fsee iBlanc etaI1[200l. zero point un- 
certainties (compare iRahman et al.l l2010t ) or a correct 
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Fig. 2. — ''"Dcp ^ function of Sh2- Gray points indicate indi- 
vidual measurements. Red points show the running median. Error 
bars indicate the la scatter in each bin. The dashed line shows the 
median r^^^ ~ 2.35 Gyr. r^^p shows little or no systematic varia- 
tion with H2 surface density over the range Sh2 ~ 5-100 M0 pc~^. 

treatment of the uncertainties associated with physical 
parameter estimation. Bearing these caveats in mind, a 
rough parameterization may still be useful to the reader. 
If we apply a simple linear regression in log space and 

fi10 the function Esfr = Ax (Sh2/10Mo pc'^)^ to 
the binned kpc data (red points in the lower right panel 
of Figure [U, we find A w 4.4 x 10"^ Mq yr"^ kpc"^ 
and N « 1.0. This is not rigorous: we have treated the 
observable I]h2 as an independent variable and we dis- 
carded information in the process of binning. However 
the fit does reasonably bisect the data. We find simi- 
lar results fitting the individual measurements where we 
are complete with N varying by ±0.2 and A varying by 
~ 30%, depending mainly on how the fit is constructed. 

The results of this fitting can be distilled to what is 
immediately apparent from the plot: a characteristic 
Tpgp ~ 2.3 Gyr and a power law index close to unity, 
so that the data extend parallel to the dashed lines of 
fixed Tp^p in Figure [T] The global index close to unity 
implies that the ratio of Eh2 to Ssfr does not change 
much as a function of Sh2 across our data. We quan- 
tify this by comparing r^^p to I]h2 where we are com- 
plete (Sh2 > 5 Mq pc~^). Figure [2] plots the individual 
measurements along with a running median and scatter; 
both show little or no systematic variation of r^^p as a 
function of SH2 across the range studied. The rank cor- 
relation coefficient relating r^^p = I]h2/Ssfr to Sh2 is 
r — 0.09 ± 0.01 in our kpc data, i.e., the two quantities 
are only very weakly correlated. 

These re s ults e xtend those found by BOS and 
iLerov et"al1 (|2008[ ). who also found a roughly con- 
stant ratio Sh2/5]sfr for a smaller, less diverse sample. 
Based on detailed studies of Local Group galaxies (e.g . 
Blitz et al.l[2007l iBolatto et al.l[200l iBigiel et al.]l2010al : 
Fukui fc Kawamurall2010[ ). thev speculated that the ap- 

We normalize the fit at Sh2 = 10 Mq pc'^ following BOS. 



proximately linear I]sfr-Sh2 relation arises because star 
formation in disk galaxies takes place in a relatively uni- 
form population of GMCs. Given typical GMC masses 
of ^ 10^-10^ Mq and sizes of ~ 50 pc, each of our res- 
olution elements likely averages over at least a few — 
and often many — GMCs. Thus, in this scenario the re- 
lationship between Eh2 and Ssfr reduces to a counting 
exercise: SH2 corresponds to a different number of GMCs 
inside different resolution elements, rather than to chang- 
ing physical conditions in the molecular gas. This also 
naturally explains the weak dependence of our results 
on spatial scale, which merely determines the number 
of GMCs per resolution element but leaves the average 
fixed constant r^^p intact (compare BOS for a detailed 
discussion). 

This scenario does not contradict earlier results find- 
ing that TQgp depends on I]h2 or Mh2- the strongest 

measurements of variable Tp^p come from LIRGs and 
ULIRGs fe.g.. IKennicut^ H^M IGao fc Soloim^ [2004). 
systems with H2 surface densities significantly exceed- 
ing those studied here and where the assumption of a 
uniform GMC population likely breaks down. Depar- 
tures are also expected on scales of individual molecular 
clouds, where only a small fraction of the molecul ar gas 
actively forms stars fe.g.. iHeiderman et al.l [2"010( ). We 
will show in the next section, however, that for normal 
disk galaxies and on scales greater than a few 100 pc 
our results agree remarkably well with previous measure- 
ments of Tpcp- 

4. COMPARISON TO LITERATURE DATA 

As described in Section[Tl many studies have examined 
the relationship between molecular gas and star forma- 
tion in nearby disk galaxies over the last decade. The 
emphasis on power law fits has somewhat obscured the 
basic question of whether these data fundamentally agree 
or disagree regarding which part of Eh2 — Ssfr space is 
occupied by local disk galaxies. To address this point. 
Figure [3] shows our binned data (big points with error 
bars) along with a wide compilation of recent measure- 
ments. 

We adjust all literature measurements to share our 
adopted Xqo and stellar IMF (Kroupa), but otherwise 
leave the data unchanged. These points therefore reflect 
a wide range of star formation tracers, sampling schemes 
and physical scales. 

We plot averages over whole galaxies as triangles. 
These include 57 normal spiral ga laxies (gre en.) and 
15 starbu r st ga laxies (blue) from iKennicut^ ([l99S). 
lKennicut3 ()1998[ ) estimates Ssfr from Ha for normal 
spirals and IR emission for sta rbursts. We al so show 236 
pointings towards spirals from iMiirgia et aP ([2002, pur- 
ple) and towards 80 small nearby spirals and dwarfs from 
ILerov et al.] (|2005L redJEJ. Both data sets have 50" 
resolution, corresponding to ^ 1-4 kpc, and use 1.4 GHz 
radio continuum (RC) emission to estimate Ssfr- 

Filled diamonds indicate radial profile measur e ments . 
Data for 7 nearby spirals from iWong fc Blitz! (I2002f) 
are shown in red, those for M51 from iSchuster et al] 

These studies comp i le measurements from |Young_et_^ 
IIT995I) , ITavlor et all ifTOOSf i. lElfhag et all pMB ). and IBoker eTal 
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Fig. 3. — SsFR versus Sh2 for a compilation of literature measurements and our binned data. In the left panel we label individual 
studies, which employ a wide range of star formation tracers, sampling schemes and physical scales. The black points indicate the running 
medians for our 1 kpc data from Figure [T] In the right panel we treat all literature measurements equally (gray points) and construct a 
running median (red points) in the same way that we binned our data (black points) in Figure [l] Both panels show excellent agreement 
between our measurements and the literature data and suggest an emerging consensus on the basic, approximately linear, SsFR ~ 2h2 
relation in nearby disk galaxies. 



(120071) in green and tho se for N GC 6946 from 
iCrosthwaite fc Turned (|2007[ ) in purple. iWong fc Blitz 
( 2002f ) derive Esfr from Ha ei nission. iSchuster et al. 
(2007) from RC emission, and ICrosthwaite fc Turner 
11007) from FIR emission. 

Small points represent aperture data. Blue 
points show 520pc-siz ed aperture measurements from 
iKennicutt et al.l ()2007[ ) of star forming regions in the spi- 
ral arms of NGC 5194 (M51). They infer Esfr from a 
combination of Ha and 24 / zm emission. Green points 
show 500 pc apertures from iRahman et al.l (|201Cl( ). who 
sample mainly the spiral arms of NGC 4254. The points 
shown here reflect Ssfr as derived from FUV and 24 fim 
emission. Red points indicate 170 pc apertures cover- 
ing the central 4. 1 x 4.1 kpc^ of NGC 5194 (M51) from 
iBlanc et a l. (2009) . They infer Ssfr from extinction cor- 
rected Ha emission using integral field unit observations. 
The left panel of FigureEllabels these various studies and 
overplots our data. 

Figure [3] shows that these measurements sweep out a 
distinct part of I]sfr~Sh2 space. Most data scatter be- 
tween Tpep = 10^ and 10^" yr and our measurements lie 
near the center of the distribution. The right panel in 
Figure [3] shows this most clearly: we take the simplis- 
tic approach of treating all of the literature data equally 
(shown as gray points) and construct the same running 
median that we use on our own data. The literature 
average (red points) agrees strikingly well with our mea- 
surements (black points). This implies that our results 
are robust with respect to the choice of tracers or ex- 
perimental setup. The literature sample as a whole also 
suggests that r^^p w 2.3 Gyr in nearby disks and that 

Tp^p is a fairly weak function of Sh2- 
5. SUMMARY 



Using new IRAM 30m CO J = 2 1 maps from 
the HERACLES survey, we determine the relation be- 
tween H2 surface density, Sh2, and SFR surface density, 
SgFR, in 30 nearby disk galaxies. This significantly ex- 
tends the number of galaxies (by more than a factor of 
four) an d the range o f gala xy properties probed com- 
pared to [BiiiiriraD (|200l) . We present our mam re- 
sults for a common physical resolution of 1 kpc. We find 
a remarkably constant molecular gas consumption time 
Tp^p ?a 2.35 Gyr (including helium) with a Icr scatter of 

0.24 dex (w 75%) and little dependence of r^^p on Eh2 
over the range I;h2 ~ 5-100 Mq pc~^. 

This ext e nds and reinforces the conclusions of 

iBigiel et all (|2008D and lLerov et al.l (pOOSl) that the star 
formation rate per unit H2 in the disks of massive star- 
forming galaxies is, to first order, constant. We interpret 
this as strong, yet indirect, evidence that the disks of 
nearby spiral galaxies are populated by GMCs forming 
stars in a relatively uniform manner. We caution that 
these results are specific to disk galaxies and scales on 
which we average over many GMCs — they may be ex- 
pected to break down at very high surface densities and 
on small scales. Taken as a whole, a broad compilation 
of literature data on disk galaxies from the last decade 
yields impressively similar results. 



We thank the GALEX NGS, SINGS, and LVL teams 
for making their outstanding datasets available. We 
thank Karl Gordon for the kernel used on the MIPS 
24/im data and Nurur Rahman for sharing his data. We 
thank the staff of the IRAM 30m telescope for their as- 
sistance carrying out the survey. F.B., A.K.L., and F.W. 
gratefully acknowledge the Aspen Center for Physics, 
where part of this work was carried out. Support for 



6 



F. Bigiel et al. 



A.K.L. was provided by NASA through Hubble Fehow- 
ship grant HST-HF-51258.01-A awarded by the Space 
Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, 



Inc., for NASA, under contract NAS 5-26555. The work 
of W.J.G.d.B. is based upon research supported by the 
South African Research Chairs Initiative of the Depart- 
ment of Science and Technology and National Research 
Foundation. 



REFERENCES 



Bigiel, F., Leroy, A., Walter, F., Brinks, E., de Blok, W. J. G., 
Madore, B., & Thornley, M. D. 2008, AJ, 136, 2846 

Bigiel, F., Bolatto, A., Leroy, A., Blitz, L., Walter, F., 

Rosolowsky, E., Lopez, L., & Plambeck, R. 2010a, ApJ, 725, 
1159 

Bigiel, F., Leroy, A., Walter, F., Blitz, L., Brinks, E., de Blok, 

W. J. G., & Madore, B. 2010b, AJ, 140, 1194 
Blanc, G. A., Heiderman, A., Gebhardt, K., Evans, N. J., & 

Adams, J. 2009, ApJ, 704, 842 
Blitz, L. 1993, in Protostars and Planets III, ed. E. H. Levy & 

J. I. Lunine, 125-161 
Blitz, L., Fukui, Y., Kawamura, A., Leroy, A., Mizuno, N., & 

Rosolowsky, E. 2007, in Protostars and Planets V, ed. 

B. Reipurth, D. Jewitt, & K. Keil, 81-96 
Boker, T., Lisenfeld, U., & Schinnerer, E. 2003, A&A, 406, 87 
Bolatto, A. D., Leroy, A. K., Rosolowsky, E., Walter, F., & Blitz, 

L. 2008, ApJ, 686, 948 
Calzetti, D., et al. 2007, ApJ, 666, 870 
Calzetti, D., et al. 2010, ApJ, 714, 1256 
Crosthwaite, L. P., & Turner, J. L. 2007, AJ, 134, 1827 
Daddi, E., et al. 2010, ApJ, 714, L118 
Dale, D. A., et al. 2007, ApJ, 655, 863 
Dale, D. A., et al. 2009, ApJ, 703, 517 

Elfliag, T., Booth, R. S., Hoeglund, B., Johansson, L. E. B., & 

Sandqvist, A. 1996, A&AS, 115, 439 
Fukui, Y., & Kawamura, A. 2010, ARA&A, 48, 547 
Gao, Y., & Solomon, P. M. 2004, ApJ, 606, 271 
Genzel, R., et al. 2010, MNRAS, 407, 2091 
Gil de Paz, A., et al. 2007, ApJS, 173, 185 

Heiderman, A., Evans, N. J., II, Allen, L. E., Huard, T., & Heyer, 

M. 2010, ApJ, 723, 1019 
Heyer, M. H., CorbeUi, E., Schneider, S. E., & Young, J. S. 2004, 

ApJ, 602, 723 
Kennicutt, Jr., R. C. 1998, ApJ, 498, 541 



Kennicutt, Jr., R. C., et al. 2003, PASP, 115, 928 
Kennicutt, Jr., R. C., et al. 2007, ApJ, 671, 333 
Leitherer, C., et al. 1999, ApJS, 123, 3 

Leroy, A., Bolatto, A. D., Simon, J. D., & Blitz, L. 2005, ApJ, 
625, 763 

Leroy, A. K., Walter, F., Brinks, E., Bigiel, F., de Blok, W. J. G., 
Madore, B., & Thornley, M. D. 2008, AJ, 136, 2782 

Leroy A. K., et al. 2009, ApJ, 702, 352 

Marble, A. R., et al. 2010, ApJ, 715, 506 

McKee, C. F., & Ostriker, E. C. 2007, ARA&A, 45, 565 

Moustakas, J., Kennicutt, R. C, Jr., Tremonti, C. A., Dale, 
D. A., Smith, J.-D. T., & Calzetti, D. 2010, ApJS, 190, 233 

Murgia, M., Crapsi, A., MoscadeUi, L., & Gregorini, L. 2002, 
A&A, 385, 412 

Rahman, N., et al. 2010, ApJ, in press. tajXiv:1009.3272i 

Rownd, B. K., & Young, J. S. 1999, AJ, 118, 670 

Salim, S., et al. 2007, ApJS, 173, 267 

Schuster, K. F., Kramer, C, Hitschfeld, M., Garcia-Burillo, S., & 

Mookerjea, B. 2007, A&A, 461, 143 
Taylor, C. L., Kobulnicky, H. A., & Skillman, E. D. 1998, AJ, 

116, 2746 

Verley, S., CorbeUi, E., Giovanardi, C, & Hunt, L. K. 2010, 
A&A, 510, A64 

Walter, F., Brinks, E., de Blok, W. J. G., Bigiel, F., Kennicutt, 

R. C, Thornley M. D., & Leroy, A. 2008, AJ, 136, 2563 
Warren, B. E., et al. 2010, ApJ, 714, 571 
Wilson, C. D., et al. 2009, ApJ, 693, 1736 
Wong, T., & Blitz, L. 2002, ApJ, 569, 157 
Young, J. S., et al. 1995, ApJS, 98, 219 

Young, J. S., Allen, L., Kenney, J. D. P., Lesser, A., & Rownd, B. 
1996, AJ, 112, 1903 



